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The efficiency of star formation in clustered and 
distributed regions 
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ABSTRACT 

We investigate the formation of both clustered and distributed populations of young 
stars in a single molecular cloud. We present a numerical simulation of a 1O 4 M0 
elongated, turbulent, molecular cloud and the formation of over 2500 stars. The stars 
form both in stellar clusters and in a distributed mode which is determined by the 
local gravitational binding of the cloud. A density gradient along the major axis of 
the cloud produces bound regions that form stellar clusters and unbound regions that 
form a more distributed population. The initial mass function also depends on the 
local gravitational binding of the cloud with bound regions forming full IMFs whereas 
in the unbound, distributed regions the stellar masses cluster around the local Jeans 
mass and lack both the high-mass and the low-mass stars. The overall efficiency of star 
formation is rs 15% in the cloud when the calculation is terminated, but varies from less 
than 1% in the the regions of distributed star formation to « 40% in regions containing 
large stellar clusters. Considering that large scale surveys are likely to catch clouds at 
all evolutionary stages, estimates of the (time-averaged) star formation efficiency for 
the giant molecular cloud reported here is only « 4%. This would lead to the erroneous 
conclusion of slow star formation when in fact it is occurring on a dynamical timescale. 

Key words: stars: formation - stars: luminosity function, mass function - globular 
clusters and associations: general. 



1 INTRODUCTION 

The ability to conduct wide-area surveys of molecular clouds 
has shown that most stars form in c lusters containing 
some hundreds to thousands of star s l|Lada et alj Il99lt 
IClarke et ai1l2000l ; lLada fc Ladal 120051 '). At the same time, 
mid-infrared surveys such as Spitzer have shown that sig- 
nificant numbers of stars form in a more distributed mode 
il Allen et al.ll2007l ; iGutermuth et alj|2008l . 120091 ; lEvans et all 
2009). The reason why such different modes of star forma- 
tion exist, and in the same cloud (e.g. Orion A) is unclear. 

There has also been considerab le interest as to wh y star 
formation appears to be inefficient (|Evans et al.l[200 9f), with 
only a few percent of a molecular cloud's mass being turned 
into stars per free-fall ti me. This could imply th at star for- 
mation is a slow process IjKrumholz fc Tanir2007) or that it 
is an inherently inefficient process, but proceeds on the local 
dynamical timescale. In the latter case the efficiency must 
increase on small scales where bound clusters are formed. 
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For example, the Orion Nebula Cluster has a median age 
of « 10 6 years and a dynamic al time of « 3 x 10 5 years 
l|Hillenbrand fc Hartmannlll998l ). Given an overall star for- 
mation efficiency of w 50% this implies a star formation 
per free-fall time of 15%. Considering that the initial pre- 
cluster cloud is likely to have been at least a factor of 2 
larger (|Bonnell et al . 2003), this implies an efficiency of star 
formation per initial free-fall time of close to 50 %. 

To date, numerical simulations have generally cho- 
sen spherically symmetric or period boxes initial condi- 
tions of gravitationally bound cl ouds which collapse and 
fragment to form stellar clusters llKlessen. Burkert fc Bate! 
1 19981 : iBate, Bonnell fc Bromml 120031 : iBatd I2009D . Cluster 
formation proceeds through hierarchical fragmentation and 
production of a somewhat distributed population which 
undergoes a hierachical merger process from small sub- 
clusters to one final c l uster containing most of the stars 
llBonnell. Bate fc Vine! 120031 : IBatd 120091 : iFederrath et~ai1 
|2010D . One simple possibility is that if star formation occurs 
in regions of molecular clouds that are globally unbound, 
then there is no reason for the stars that form from the frag- 
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meriting population to fall together to form the large stellar 
cluster. Recent work evaluating the boundness of molecu- 
lar clouds show that their masses are typically five times 
smaller than that to be virialised, implying that much of the 
present day star formation is occurring in unbound molec- 
ular clouds(?). Here we demonstrate that the outcome of a 
distributed or clustered population can depend on whether 
the region is, or is not, globally bound. 

Gravitationally unbound clouds have been explored in 
a series of studies to inv estigate how this relates to the effi- 
ciency of star formation dClark fc Bonnelj[2003 ;l Clark et al.l 
120051 : IClark. Bonnell fc Klessenl l2008t ). Low star formation 
efficiencies are commonly taken to imply that star forma- 
tion is slow and that molecular clouds are long-lived entities, 
supported by some internal mechanism and lasting for sev- 
eral tens of dynamical times. In contrast, unbound clouds 
can also produce low star formation efficiencies on dynam- 
ical timescales due to the fact that only a fraction of the 
cloud becomes gravitationally bound due to the turbulence 
and undergoes gravitational collapse and star formation. 

In this paper, we explore the importance of the local 
gravitational binding in one cloud and show that a single 
cloud can produce both a distributed and a clustered popu- 
lation, and a range of star formation efficiencies, depending 
on the local gravitational binding. 



2 CALCULATIONS 

The results presented here are based on a large-scale 
Smoothed Particle Hydrodynamics (SPH) simulation of a 
cylindrical 10 4 Mq molecular cloud 10 pc in length and 

3 pc in cylindrical diameter. We have chosen an elon- 
gated cloud rather than the more standard spherical cloud 
as most molecular clouds are non-sperhical and commonly 
elongated (e.g. Orion A). Such a geometry can also pro- 
duce additional structure due to gravitational focussing 
l|Hartmann fc B urkcrt 20q3). This also allows for the phys- 
ical properties to be varied along the cloud in a straight- 
forward manner. The cloud has a linear density gradient 
along its major axis with maximum/minimum values, at 
each end of the cylinder, 33 percent high/lower than the 
average gas density of 1.35 x 10~ 20 g cm -3 . The gas has in- 
ternal turbulence following a Larson-type P(k) ~ k~ 4 power 
law throughout the cloud and is normalised such that the 
total kinetic energy balances the total gravitational energy 
in the cloud. This corresponds to a full cloud (10 pc) 3-D 
velocity dispersion of order 4.5 km s . The density gradi- 
ent applied then results in one end of the cloud being over 
bound (still super virial) while the other end of the cloud is 
unbound. 

The cloud is populated with 15.5 million SPH parti- 
cles on two levels, providing high resolution in regions of 
interest. We initially performed a lower resolution run with 
5 million SPH p articles producing an average mass reso- 
lution of 0.15Mm iBate fc Burkertl (|l997T l. Upon completion 
of this low resolution simulation, we used three criteria to 
identify the regions that required higher resolution. This 
included the particles which formed sinks, and those that 
were accreted onto sinks. It also included particles which at- 
tained sufficiently high density such that their local Jeans 
mass was no longer resolved in the low-resolution run. All 



of these particles were identified and from the initial condi- 
tions of the low resolution run, they were split into 9 par- 
ticles each to create the initial conditions for the high reso- 
lution simulations. This particle splitting was performed on 
the initial conditions to ensure that the physical quantities 
of mass, momentum, energy and the energy spectrum were 
preserved. Note that the particle splitting does not intro- 
duce finer structure in the turbulent energy spectrum. This 
produced a mass resolution for the regions involved in star 
formation of 0.0167Mq , sufficient to resolve the formation of 
higher-mass brown dwarfs, equivalent to a total number of 
4.5 x 10 7 SPH particles. The equation of state (below) was 
specified in order to ensure that the Jeans mass in the higher 
resolution run did not descend below this mass resolution. 

Particle splitting results in a marked increase in 
resolution without unmanageable computational costs 
l|Kitsionas fc Whitworthll2002ll2007l '). Note however some of 
the unsplit particles, which in the low resolution run neither 
exceeded their Jeans mass limit nor became involved in the 
star formation, did get accreted by the additional stars in 
the high resolution run. This is to be expected as there are 
now additional locations of of star formation not present in 
the low resolution run and these additional sinks will neces- 
sarily accrete unsplit particles. 

Th e simula t ion fo llows a modified Larson-type equation 
of state iLarsonl l|2005h comprised of three barotropic equa- 
tions of state 

P = kp 1 (1) 
where 
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and p x = 5.5 x 10" 19 g cm~ 3 ,p 2 = 5.5 x 10" 15 g cm~ 3 ,p 3 = 
2 x 10 _13 g cm" 3 . 

The initial cooling part of the equation of state mimics 
the effects of line cooling and ensures that the Jeans mass at 
the point of fragmentation is appropriate for charact eristic 
stellar mass (Ijappsen et alj|200a bonnel! et alj|2006j ). The 
7 = 1.0 approximates the effect of dust cooling while the 
7 = 1.4 mimics the effects of an optically thick (to IR radi- 
ation) core, although its location at p — 5.5 x 10~ 15 g cm -3 , 
at lower densities than is typical, is in order to ensure that 
the Jeans mass is always fully resolved and that a single self- 
gravitating fragment is turned into a sink particle. A higher 
critical density for this optically-thick phase where heating 
occurs would likely result in an increase in the numbers of 
low mass objects formed. The physical processes described 
would be unchanged. The final isothermal phase of the equa- 
tion of state is simply in order to allow sink-particle forma- 
tion to occur, which requires a subvirial collapsing fragment. 
The initial conditions of the cloud contain 891 thermal Jeans 
masses (Mj oans ~ 11M©) such that if the cloud were isother- 
mal, we would expect of order 900 fragments to form. 

Star formation in the cl oud is modelled through th e 
introduction of sink-particles l|Bate. Bonnell fc Priceill995l ). 
Sink-particles formation is allowed once the gas density of 
a collapsing fragment reaches p > 6.8 x 10 -14 g cm -3 
although the equation of state ensures that this requires 
p ^ 2.x 10~ 13 g cm -3 . The neighbouring SPH particles need 
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Figure 1. The evolution of our model GMC as it evolves to form both a distributed (bottom) and clustered (top) population of stars. 
The cloud is initially globally bound but a density gradient along the major axis of the cloud makes the lower region unbound while the 
top region is gravitationally bound. The cloud is shown at 0.365, 0.544, 0.727 0.961 tg (tff m 6.6 X 10 5 ) . Each panel shows the cloud in 
a lOxlOpc region. The gas column densities are plotted on a logarithmic scale from 0.01 (black) to 100 (white) g cm -2 . 



be within a radius of 1. x 10~ 3 pc and that fragment must 
be subvirial and collapsing. Once created, the sinks accrete 
bound gas within 1. x 10 -3 pc and all gas that comes within 
2. x 10 -4 pc. The sinks have their mutual gravitational inter- 
actions smoothed to 2. x 10~ 4 pc or 40 au. No interactions 
including binary or disc disruptions can occur within this 
radius. 

We assume a 100 % efficiency of star formation within 
our sink particles. This likely overestimates the efficiency 
that would result were feedback from massive stars included. 
It is worth noting that our gas densities and core sizes are 
similar to the continuum surveys (Andre ..) that would re- 
quire a 100 % conversion in order to obtain a mapping from 
the core mass function to the stellar IMF. Furthermore, Pre- 



vious simulations including ionisation and winds l|Dale et al.l 
l2005l ; |Pale fc BonnelfeOOSI ) do not find a large change in the 
resultant masses or mass spectra. 



3 STAR FORMATION AND THE 
DEVELOPING IMF 

The simulation was followed for 1.02 free-fall times or ~ 
6.6 x 10 5 years and ~ 3.9 x 10 years after the first stars 
formed (see Fig.[TJ. During this time, 2542 stars were formed 
with masses between 0.017 and 30 Mq . The majority of these 
stars form in the upper gravitationally bound part of the 
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Figure 2. The developing initial mass function of our star forming GC is plotted at various times given in units of the clouds free-fall 
time, tg w 6.6 x 10 5 . 



cloud while some 7 per cent form in the lower, gravitationally 
unbound regions. 

Figure [2] shows the developing initial mass function 
during the star formation process. The stars form with 
masses comparable to the Jeans mass of the local gas. 
These initial masses are initially of the order of several 
tenths of a Solar mass, while lower mass fragments form 
stars later in the evolution due to the compression of gas 
to h igher densities as it falls into existing stellar clus- 
ters (|Bate. Bonnell fc Bromml [200^ ; iBonnell. Clark fc Bate! 
2008). Low and intermediate mass stars located in the cen- 
tre of forming clusters continue to acc reted from the infalling 
gas and become high-mass stars (e . g. iBonnell. Vine fc Bate! 
12004 ISmith. Longmore fc Bonnelll I2009T ). This produces a 
mass function that resembles the stellar IMF at all points 
during the evolution with a continuous source of low mass 
stars forming with a decreasing subset of these accreting 
to ever higher masses. The hig h-mass end of the IMF i s 
somewhat natter than Salpeter l|Maschberger et al.ll201ol ). 
This leaves room for the additional physics of feedback from 
massive stars and the expected decrease in efficiency of mas- 



sive star formation. Note that by the massive stars attain 
their high-mass status t hrough ongoing accr etion over rela- 
tively long time-periods (|Bonnell et al.ll2004h such that their 
feedback could only affect the cloud after much of the star 
formation has occurred. 



The cloud produces a variety of outcomes in terms of 
the distribution of stellar masses, clustered and distributed 
modes of star formation as well as the efficiency of the star 
formation process. These all depend largely on the initial 
conditions of the cloud and in particular to how gravitation- 
ally bound the cloud is locally. The density gradient that is 
imposed along the major axis, in conjunction with a con- 
stant specific kinetic energy of the gas, results in a local 
variation of the gravitational binding. Measured in terms of 
the critical mass per unit length to be bound, this variation 
extends from M/L = 0.6 (unbound) to M/L = 1.4 (bound 
with the cloud overall having a M/L = 1. 
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Figure 3. The formation of stellar clusters occurs as infalling gas filaments fragment to form small-N clusters, these clusters grow through 
the accretion of gas and stars and the merger of subclusters. The cloud is shown at 0.58, 0.73, 0.87 and 1.02 % (% =fB 6.6 X 10 5 ). Each 
panel shows the cloud in a 2.0x2.0pc region. The gas column densities are plotted on a logarithmic scale from 0.01 (black) to 100 (white) 
g cm -2 . 



4 CLUSTERING 



The evolution produced a number of high density clusters as 
well as a distributed population of stars. The clusters form 
predominantly in the (upper) bound regions of the cloud. 
The clusters form through the fragmentation of local over- 
dense filamentary structures that arise due to the turbu- 
lence, especially where such filaments intersect. Stars fall 
into local potential wells and form small-N clusters which 
quickly grow by accreting other stars (and gas) that flow 
along the filaments into the cluster potential. The merger of 
clusters also contributes to the growth of a stellar cluster. 
Figure [3] shows an example of this process whereby accre- 



tion of gas and stars occurs along filaments flowing into the 
cluster. 

An important result from this work is that the cluster- 
ing depends strongly on the local gravitational binding of 
the gas prior to star formation. Figure [4] shows the resulting 
stellar densities as a function of how bound the cloud was 
initially in terms of the critical mass per unit length, M/L, 
to be bound. Two measures of the stellar densities are plot- 
ted. The blue crosses show the median stellar density de- 
termined by the volume needed to contain the ten nearest 
stellar neighbours. The black triangles show the density of 
stars contained in a fixed volume of size 0.5pc. The density 
determined by the first method is significantly higher as it 
typically is based on much smaller volumes. In both cases, 
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Figure 4. The stellar density is plotted as a function of the local 
gravitational binding measure in terms of the critical mass to 
length measure for a cylinder to be gravitationally bound. The 
higher (blue) points show the median density measured from the 
ten nearest neighbours to each "star" while the lower value is 
measure over the local volume of the cloud. We see that low stellar 
densities correspond to unbound regions of the cloud whereas high 
stellar densities result from bound regions. 



the stellar density is low in regions that were initially un- 
bound and is much higher in the bound parts of the cloud 
with M/L > 1. 

This result is understandable in that clusters are (at 
least temporarally) bound objects and their formation re- 
quires that the pre-star formation gas is also bound. In lo- 
cally unbound regions, it is still possible to form small stellar 
clusters in regions where turbulent compression and shocks 
result in a locally bound region. This process is more effi- 
cient in regions that are globally bound. Larger scale regions 
containing many subsystems are bound even before any tur- 
bulent support is dissipated. The smaller systems that form 
locally c an then hierarchical l y merge to form large stellar 
clusters l|Bonnell et al.l 120031 : iBatd l2009h . Residual gas in 
these bound regions then falls into the gravitational po- 
tential of the cluster to be competitively accreted by the 
growi ng massive stars located in the bottom of t he potential 
well (|Bonnell et al.ll2004l ;l Bonnell fc Batd fcoOG). The mas- 
sive stars are thus located in the stellar clusters. 

The majority of the stars and brown dwarfs formed are 
in high-density regions or have been ejected from s tellar clus- 
ters thmtigli_jrr^actions (|Bate et al.ll20oj .r2003). As noted 
in (|Bonnell et al.l [20081 ) . the brown dwarfs predominantly 
form in stellar clusters due to the compression of the gas to 
high local densities as it falls into the gravitational poten- 
tial. The high-mass stars are also predominant ly formed in 
clusters l|Bonnell et al.1 12004 ISmith et al.ll2009l ). This leads 
to a potentially observable difference in the stellar IMFs of 
distributed and clustered star formation. Figure [5] shows the 
final IMF for the overall population and also for distributed 
and clustered populations defined as those with a stellar 




Figure 5. The initial mass functions at the end of the simulation 
are plotted for the total (top, black line) population and for those 
stars formed in high-density regions (over 100 stars pc -3 , red 
line) and low density regions (lower, blue line). We see that the 
low density regions do not form either high or low mass stars and 
thus result in an unusual IMF. The right hand panel shows the 
same three distributions but where the stellar density is measured 
at the end of the simulation. Several stars formed in high-density 
regions have been ejected from the clusters into the field. 
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Figure 6. The cumulative mass functions are plotted for the total 
(solid black line) population and for those stars formed in high- 
density (over 100 stars pc — 3 , short-dashed red line) and low den- 
sity (long-dashed blue line) regions. The right hand panel shows 
the same three distributions but where the stellar density is mea- 
sured at the end of the simulation. We see the low-density pop- 
ulation is significantly different from that formed in the stellar 
clusters. This difference remains even after stars are lost from the 
clusters into the low-density regions. 



density lower or higher than 100 stars pc . The right-hand 
panel of figure [5] shows the corresponding populations where 
they are separated by their maximum stellar density during 
their evolution. The cumulative distributions (figure[6]) show 
that the two distributions are statistically different and in- 
consistent (at the 1 x 10 -13 level) with being drawn from 
the same population. The distributed population has a sig- 
nificantly higher median stellar mass and a pronounced lack 
of low-mass objects. This result helps explain the seemingly 
anomolous IMF in Taurus which appears to have a lack of 
brow n dwarfs and h igh-mass stars in a distributed popula- 
tion l|Luhmanll2003 'l. 



5 EFFICIENCY OF STAR FORMATION 

One of the central questions we wish to address in this 
paper is the relationship between the nature and effi- 
ciency of the star formation process. Previou s studies (i.e. 
IClark fc Bonnellll2004 rGlark et al.ll2005l , 12008ft showed that 
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Figure 7. The star formation efficiency is plotted as a function 
of the critical mass to length ratio for a cylinder to be gravita- 
tionally bound. High star formation efficiencies result from very 
bound regions whereas very low star formation efficiencies result 
in unbound regions. A mixture of such regions can thus result 
in the formation of high SFE clustered regions and low SFE dis- 
tributed populations. 



unbound clouds resulted in inefficient star formation. Fur- 
thermore, the efficiency reduces dramatically the further the 
clouds are from being bound. Turbulent compression and 
shocks results in some star formation in these clouds but it 
is localised and much of the cloud escapes without entering 
the star formation process. 

In the present study, we have one cloud that has regions 
which are bound and regions which are unbound with a spa- 
tially varying M/L from 0.6 to 1.4. This results in a range 
in local star formation efficiencies from 0.006 to 0.4. Fig- 
ure [7] plots the local star formation efficicency as a function 
of the local binding of the cloud in terms of the critical M/L 
for the cloud to be globally bound. We see that after 1.02 
free-fall times or « 6.6 x 10 5 years (and « 3.9 x 10 5 years 
after the first stars formed) the local efficiency of star for- 
mation is strongly dependent on local binding of the cloud. 
The bound regions have efficiencies varying from 10% where 
the cloud is just bound to 20% for a M/L = 1.2 to a peak 
value of 40% near where the cloud is maximally bound. On 
the unbound side, the efficiency quickly drops below 10%, 
reaching values as low as a few % for M/L < 0.8. We can 
thus conclude that small changes in the local binding of the 
cloud result in vastly different outcomes in terms of the star 
formation efficiency. 

It is worth noting that there is a strong correlation be- 
tween the local efficiency of star formation and the formation 
of stellar clusters. Regions with relatively high efficiencies of 
^ 10% corresponds to regions which are bound and thus 
form stellar clusters. In contrast the unbound regions form 
a relatively distributed, low stellar density population and 
does so at very low efficiencies. This is in agreement with ob- 
servations where clustered regions are found to have higher 



star forming efficiencies, whereas distributed regions such 
as Taurus have low star formation efficiencies. These dif- 
ferences in the star formation efficiencies reported here are 
not simply due to delays in star formation in the unbound 
parts of the cloud as most of the mass is actually leaving 
the cloud and can not partake in the star formation process 
l|Clark et alj|2008h . 



6 OBSERVABLE DETERMINATION OF THE 
STAR FORMATION EFFICIENCY 

The evolution of the star formation efficiency is shown in 
figure [8] from where the first stars form at t m 0.5tft to the 
end of the simulation at t « l.Otg. The overall star formation 
efficiency as measured at the end of the simulation is ~ 15 %. 
this global val ue is an upper limit as no feedback effects are 
included (e.g. iDale fc Bonnellll2pjosh . Magnetic fields could 
also act to reduce this number further l|Price fc Batell2008l . 
2009). For example, if feedback acted to destroy the cloud 
at 0.7tft , then the final star formation efficiency would be of 
order 5%. 

Such global star formation efficiencies are commonly in- 
voked to discriminate between slow, and fast star formation. 
Slow star formation invokes some supporting mechanism to 
prolong the lifeti me of molecular clouds to many tens of 
dynamical times (|Krumholz fc Tanll2007t) whereas fast star 
formation is ex pected to occur o n timescales of several dy- 
namical times pimegreenlbOOOl ). Comparing such ideas to 
observations of large scale star formation rates can be prob- 
lematic as we do not know if the measured star formation 
efficiencies are final, intermediate or even pre-star formation 
values. Figure [5] also shows the time averaged star forma- 
tion efficiency at any given time during the cloud's evolu- 
tion. This takes into account that when measuring global 
volume averaged values, we are just as likely to observe any 
particular cloud at any point in its evolution. The time av- 
eraged star formation efficiency is significantly lower than 
the instantaneous value throughout its evolution, due to the 
amount of time the cloud spends in its pre-star formation 
stage. This neglects the time take for the cloud to form which 
should be of order or greater than the dynamical time of the 
cloud. Nevertheless, the final time-averaged star formation 
efficiency is « 3 — 4%. So, in a volume where a mixture of 
clouds of different evolutionary stages are present, an ob- 
server would estimate star formation efficiency of several 
percent and thus conclude that this was slow star forma- 
tion, when in reality star formation is locally proceeding 
on a dynamical timescale. All of the above neglects the ef- 
fects of magnetic fields and feedback which both act to re- 
duce the star formation rates and effici encies. While mag- 
netic fields do slow down star formation l|Price fc Bate! 2008. 
2009), feedback impedes or locally stops star form ation with- 
out appreciably changing i ts dynamical nature (|Dale et al.l 
l2005l ; lDale fc Bonnellfeoosh 



7 DISCUSSION: THE EVOLUTION OF GIANT 
MOLECULAR CLOUDS 

We have seen from the above that low star formation ef- 
ficiencies are plausibly due to the GMCs not being com- 
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Figure 8. The total (dashed) and time-averaged (solid) star for- 
mation efficiencies are plotted as a function of time. While the 
final instantaneous SFE is of order 15 %, an observer detecting a 
full population of such clouds would likely sample all evolution- 
ary stages and thus measure the much lower time-averaged star 
formation efficiency of i=a 4%. 



pletely bound by their self gravity. Even in the absence of 
effects such as magnetic fields and feedback, clouds that 
are globally, or at least in large part unbound result in low 
star formation efficiencies. Star formation requires that the 
clouds are nearly bound or at least have significant regions 
which are close to being bound, with typically gravitational 
and kinetic energies within a fa ctor of a few of each other 
l|Clark. Bonnell fc K lesscn 2008|). From this starting point, 
we can try to construct an evolutionary scenario for GMCs 
in which self-gravi ty only plays a role in the actual star 
formation process (|Pringle et al.ll200ll ; iBonnell et al]|2006l : 
iDobbs et al.ll2006h . We will neglect the effects of the mag- 
netic field but to first order its effects would be to increase 
the internal energy of the clouds such that they are further 
from being gravitationally bound. 

The formation of giant molecular clouds is uncertain 
with various theories from grav itational instabilities to c loud 
coagulation and spiral shocks l|McKee fc Os trikcr 200?]) ■ In 
reality all of these processes may play a significant role in 
GMC formation but for our purposes, we will assume that 
GMC formation oc curs primarily due to spiral shocks or 
cloud coagulations (|Dobbs et al.l r2006). In such cases, self- 
gravity need not play an important role in the formation 
process and clouds can be formed without being gravita- 
tionally bound. 

Cloud formation from external collisions or compres- 
sion implies an increasing contribution to the cloud's gravi- 
tational potential. During formation, the cloud evolves from 
a state where self-gravity is unimportant to one where self- 
gravity has a significant effect on the cloud's dynamics. Once 
the gravitational energy is within a factor of a few of the ki- 
netic energy, then star formation proceeds in local regions 
that become bound due to the cloud's internal dynamics 



l|Clark. Bonnell fc Klesser]|200Sl ). This local star formation 
will occur as parts of the cloud are still being assembled as 
the local timescale is much shorter than the overall dynam- 
ical time for the cloud or its precursor. 

Star formation would then proceed until either the local 
gas reservoir is depleted, or until the gas reservoir is removed 
by the effects of feedback from young stars. The tidal shear 
from leaving the spira l arm potential c ould also limit the 
lifetime of the clouds (|Dobbs et alj|200fj ). The majority of 
the cloud need never become gravitationally bound before 
the cloud is dispersed resulting in inefficient star formation 
process that is still occurring on a fast dynamical timescale. 



8 CONCLUSIONS 

Star formation in realistic GMCs will proceed from a variety 
of physical conditions, spanning regions that are gravitation- 
ally bound to parts or whole clouds which are gravitationally 
unbound. Star formation will occur as long as the local con- 
ditions are close to being gravitationally bound but the prop- 
erties of the young stellar population can depend strongly 
on these conditions. Regions that are bound produce bound 
stellar clusters and a stellar population that follows the full 
initial mass function from brown dwarfs to high-mass stars. 
Regions that are unbound are likely to produce a somewhat 
skewed IMF biased towards the local Jeans mass in the gas 
and with s ignificant lack of lower-mass stars, such as is seen 
in Taurus (|Luhmanll200i ). 

The star formation efficiency is also a product of the 
local physical conditions with bound regions resulting in a 
relatively high star formation efficiency of order 10 % or 
more per free-fall time. Regions that are unbound can have 
drastically reduced efficiencies of order 1 % or less per free- 
fall time. Thus clustered star formation should occur in re- 
gions of higher local star formation efficiencies that more 
distributed populations. 

Estimates of low star formation effiiciencies are equally 
consistent with fast dynamical star formation as slow qua- 
sistatic star formation provided that one relaxes the condi- 
tion that GMCs are globally bound long-lived entities. In- 
cluding the pre-star formation timeperiods where clouds are 
being assembled, global estimates of depletion timescales or 
star formation rates per free-fall time will appear to be low 
even while local regions are undergoing fast star formation 
at high efficiencies. 

Finally, realistic GMCs are likely to be constructed from 
a mix of physical conditions such that a fraction of the 
cloud is bound producing stellar clusters at high efficiencies 
whereas the majority of the cloud is unbound producing a 
more distributed population at low star formation efficien- 
cies before the cloud is unbound by feedback or alternative 
process. Such a scenario is consistent with a model where 
GMCs are not formed due to their self-gr avity but rather to 
an external p rocess such as spiral shocks l|Dobbs et al.ll2006l ; 
iDobbslbOOcf ). 
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